Coexistence of ferromagnetism and superconductivity 
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A comprehensive theory is developed that describes the coexistence ol p-wave, spin-triplet super- 
conductivity and itinerant ferromagnetism. It is shown how to use field-theoretic techniques to 
derive both conventional strong-coupling theory, and analogous gap equations for superconductivity 
induced by magnetic fluctuations. It is then shown and discussed in detail that the magnetic fluc- 
tuations are generically stronger on the ferromagnetic side of the magnetic phase boundary, which 
substantially enhances the superconducting critical temperature in the ferromagnetic phase over 
that in the paramagnetic one. The resulting phase diagram is compared with the experimental 
observations in UGe2 and ZrZn2. 
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I. INTRODUCTION 

It is well known that, in principle, the exchange of 
magnetic fluctuations between electrons can induce su- 
perconductivity in both the paramagnetic and ferromag- 
netic phases of metals.^ In general, for superconductivity 
to occur one needs both some sort of attraction between 
quasi-particles, which can be provided by magnetic flu- 
cutations, and low temperatures. Since magnetic fluc- 
tuations become large near continuous magnetic phase 
transitions, ideal candidates for this phenomenon would 
seem to be itinerant ferromagncts with a low Curie tem- 
perature. In contrast to the much more common phonon- 
exchange case, which usually leads to spin-singlet, s-wave 
superconductivity, the magnetically mediated pairing is 
believed to usually be strongest in the spin-triplet, p- 
wave channel. This type of superconductivity is very 
sensitive to nonmagnetic disorder, so that very clean sam- 
ples are also required. The combined requirements of low 
temperatures, high purity, and vicinity to a continuous 
ferromagnetic transition severely restricts the number of 
materials where ferromagnetically induced superconduc- 
tivity might be observed. Furthermore, recent theoretical 
work has shown that at low temperature the ferromag- 
netic transition in an itinerant electron system should be 
generically of first order, ^■'^ and there is some experimen- 
tal evidence for this as well.^ If this is the case, then 
there are no divergent magnetic fluctuations even at the 
transition point itself. 

Very recently, however, coexistence of ferromagnetism 
and superconductivity has indeed been observed, in 
UGe2.^ The f-electron structure of the uranium notwith- 
standing, this material is more similar to d-band metals 
than to heavy fermion systems, and ferromagnetism and 
superconductivity are believed to be caused by itinerant 
electrons in the same band. The persistence of ferromag- 



netic order within the superconducting phase has been 
ascertained by neutron scattering. Since superconduc- 
tivity in the presence of ferromagnetism is likely to be of 
spin-triplet type, magnetic-fluctuation induced pairing is 
a possible mechanism. However, two aspects of the ex- 
periments are not readily reconcilable with previous theo- 
ries of superconductivity in ferromagncts. First, the mag- 
netic transition is observed to be of first order, so that 
there are no divergent fluctuations as the transition is 
approached. Second, the superconducting state is found 
only on the ferromagnetic side of the phase boundary. A 
schematic phase diagram for this situation is shown in 
Fig. 1(a). The theoretical prediction, on the other hand, 
has been a phase diagram like the one shown in Fig. 1(b), 
with superconductivity existing, more or less symmetri- 
cally, on both sides of the magnetic phase boundary.^'^ 

Since the initial experimental observation of supercon- 
ductivity in UGe2 there has been considerable new ex- 
perimental and theoretical work. Experimentally, the 
qualitatively same phenomenon has been observed in 
ZrZn2,^ and in URhGe.^ In ZrZn2, the paramagnetic- 
to-ferromagnetic phase transition appears to be either 
continuous or only weakly first order; in URhGe the na- 
ture of the magnetic transition at very low temperatures 
is not known. In both of these materials, superconductiv- 
ity is observed far inside the ferromagnetic phase. This 
is in contrast to UGe2, where it is observed only close 
to the magnetic phase boundary, see Fig. 1(a). Theoret- 
ically, Shick and Pickett^" have included correlation ef- 
fects in a local density approximation approach and have 
concluded that magnetically mediated triplet supercon- 
ductivity can occur in UGe2. In contrast, Shimahara 
and Kohmoto^^ have argued that triplet superconduc- 
tivity can be induced by screened phonon interactions in 
ferromagnetic compounds such as UGe2. Santi et al.-'^^ 
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FIG. 1. Schematic phase diagram showing the paramag- 
netic (PM) , ferromagnetic (FM) , and superconducting phases 
(SC) in a temperature (T) - control parameter (CP) plane, 
(a) shows the qualitative PM-FM phase boundary (solid line) 
as well as the qualitative superconducting phase boundary 
line as observed in UGe2, Ref. 5 (solid line) and ZrZn2, Ref. 8 
(dashed line), respectively, and explained by the theory pre- 
sented here. In Ref. 5, hydrostatic pressure serves as CP. (b) 
shows the qualitative prediction of paramagnon theory, Ref. 
6. 

have performed ab initio calculations for ZrZn2, conclud- 
ing that magnetically mediated triplet superconductiv- 
ity can occur, with a critical superconductivity temper- 
ature as high as 1 K. Watanabe and Miyake^^ have ar- 
gued that coupled charge density waves and spin density 
waves can be used to understand the superconductivity 
in UGe2, and a Larkin-Ovchinnikov-Fulde-Ferrell state 
is also a possibility.^^ Machida and Ohmi,^^ and more 
recently Samokhin and Walker, have investigated pos- 
sible order parameters for magnetically induced super- 
conductivity, and a d-electron exchange mechanism has 
been proposed. Finally, the experiments have prompted 
investigations of the detailed role of the normal self en- 
ergy in the Eliashberg equations for magnetically induced 
superconductivity, ^^"^o 

In a previous short paper, the current authors and 
collaborators have shown that the critical temperature 
for spin-triplet, p-wave superconductivity mediated by 
spin fluctuations is generically much higher in a Heisen- 
berg ferromagnetic phase than in a paramagnetic one, 
due to the coupling of the magnons to the longitudi- 
nal magnetic susceptibility, and that this qualitatively 
explains the observed phase diagram. In the current 
paper we explain and expand on this work. There are 
two important aspects. The first, formal one, is to show 
how to use field-theoretic techniques, rather than many- 
body perturbation theory, to derive Eliashberg-like equa- 
tions for magnetic-fluctuation induced superconductivity. 
For completeness, and to motivate the more complicated 
case of magnetic-fluctuation induced superconductivity, 
we also show how these ideas can be used to derive con- 
ventional Eliashberg theory. The second one is a detailed 
explanation of our previous observation that magnetic 
fluctuations are generically stronger on the ferromagnetic 
side of the phase boundary, which in turn is shown to lead 
to a much higher critical temperature for superconductiv- 



ity in the ferromagnetic phase. This result is consistent 
with the observed phase diagrams in UGe2 and ZrZn2, 
schematically shown in Fig. 1(b). The formalism devel- 
oped here will also lay the ground for further investiga- 
tions of the superconductor-ferromagnet phase diagram. 
In particular, we will argue that as a function of temper- 
ature and magnetization, one expects phase transitions 
to other, more exotic, superconducting states inside the 
superconducting phase. 

The plan of this paper is as follows. In Section II we 
derive an effective action or field theory for superconduc- 
tivity due to magnetic fluctuations. The field theory is 
made explicit by keeping magnetic fiuctuations to one- 
loop order. The resulting superconducting field theory is 
then solved in a saddle-point approximation in Section 
III. The resulting gap equations are consistent with ear- 
lier results,^ except that the interaction is an exact lon- 
gitudinal magnetic susceptibility rather than an approx- 
imate one. The field theoretic derivation also has struc- 
tural advantages that we plan to exploit in the future, 
as is discussed in Section V. In the last part of Section 
III, explicit expressions for the magnetic susceptibility in 
either phase are given. Emphasis is on the difference be- 
tween magnetic fluctuations in the two phases. In Section 
IV the linearized gap equations derived in Section III are 
solved using a McMillan-type approximation. Phase dia- 
grams are derived showing the paramagnetic, ferromag- 
netic and superconducting phases. It is concluded that 
the critical temperature for superconductivity is much 
higher in the ferromagnetic phase than in the paramag- 
netic phase. The paper is concluded in Section V with a 
discussion of our results, of possible future applications 
of the formalism developed here, and of some open ques- 
tions. The derivation of conventional Eliashberg theory 
by field-theoretic methods is relegated to an appendix. 

II. COUPLED FIELD THEORY FOR NORMAL 

STATE, SUPERCONDUCTING, AND 
MAGNETIC FLUCTUATIONS 

A. The model 

Our starting point is an interacting electron gas model. 
In terms of Grassmann variables, the partition function 
is 

Z = J D[i;,i;] e^I^'^^l . (2.1) 

Here the functional integration is with respect to 
Grassmann- valued fields •ip and tjj, and the action S reads 

S = - f dx^ijjaix) dr V'a {x) +So + Sint ■ (2.2a) 

We use a {d + l)-vector notation, with x = (x, r), and 
J dx = Jy dx dr. x denotes position, r imaginary 



2 



time, y ^ oo the system volume, and (3 = 1/T the 
inverse temperature, a is a spin label, and we use units 
such that ?i = fee = 1. describes free electrons with 
chemical potential fi, 



^0 = y dx'^tpa{3 



2me 



M V'<T(a;) 



(2.2b) 



with TTio the electron mass. S'int denotes the two-body 
electron-electron interaction. The underlying Coulomb 
interaction becomes screened, and after renormalization 
from a microscopic to a mesoscopic scale, S'int sepa- 
rates into spin-singlet and spin-triplet particle-hole and 
particle-particle interactions, respectively.^^ The interac- 
tion in the particle-particle channel is what causes su- 
perconductivity, and we do not include it in our bare 
action. Rather, one of the main points of this paper is 
to see how this interaction is generated by magnetic fluc- 
tuations. The interactions in the particle-hole channel 
are 



with. 



SP-^ = ^ J dxn{x)n{x) 
St'"^ ^ Y / "^^^^ ■ "^^^^ 



(2.2c) 

(2.2d) 
(2.2e) 



Here Fg^t are the particle-hole spin-singlet and spin- 
triplet interaction amplitudes, respectively. 



(2.3a) 



is the electron number density field, and 
ni{x) = ^Mx)afMx) , (i = l,2,3) , (2.3b) 



are the components of the electron spin density field ng. 
o'i,2,3 denote the Pauli matrices. 

Below we will often use a Fourier representation of the 
Grassmann field. 



V'<T,n(k) 



V'cr,n(k) = 



J dx e'^'^"^-''-") i;„{x) , (2.4a) 



dx e 



(2.4b) 



Finally, to introduce the relevant order parameters, it will 
be convenient to introduce generalized Nambu spinors. 
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/ Vt(3;) 



*(a;) = -^ (V'T(a;),V'i(a;),V'T(a;),V'i(a;)) 



(2.5a) 
(2.5b) 



and corresponding composite variables as the tensor 
product of the spinors, 

^ij{x,y) = ^i{x)-^j{y) . (2.5c) 

It is also useful to define a four-vector of 4 x 4 matrices 

(70,7) = (O'S <?> Cro,(T3 (S) a-i,iTo (S) 0-2, 0-3 (g) 173 ) , (2.6) 

with (To the 2x2 unit matrix. In terms of these spinors 
and the 7, the noninteracting part of the action reads 

S - Sint = j '^^ [-^r + 70 (vV2me + n)] ^{x) , 

(2.7a) 

and the number density and spin density fields that make 
up the interacting part can be written 

n{x) = ^{x)j(,<i{x) , (2.7b) 
n:(x) = *(x)7,*(a;) , (z=l,2,3) , 

(2.7c) 

with the product between spinors defined as the matrix 
product. 



B. Superconducting and magnetic order pEirameter 
field theory 

1. Intuitive argument for magnetic-fluctuation induced 
superconductivity 

To motivate our technical considerations, we recapit- 
ulate a simple physical argument given by De Gennes.^^ 
Let us suppose that we have a local magnetization fluctu- 
ation, (5M(x), that couples to the electronic spin density. 
Since the coupling of SM. to iig is analogous to that of 
a magnetic field, the change in energy of the electron 
system is given by a Zeeman term, 

5H = -T j dx ns(x) • (5M(x) , (2.8a) 

with r a coupling strength, and ns(x) the electronic spin 
density. To first order in F, the magnetization induced at 
point X by an electronic spin density at point y is given 

by 



5Mi{x] 



Xi,(x-y)n^'(y) , (2.8b) 



where Xij is the electronic magnetic susceptibility tensor. 
This implies 



SH = -F 



' f dxdyJ2K{x)xij{^-y)ni{y) . (2.9) 
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Equation (2.9) shows that there is an effective electronic 
spin-spin interaction potential tensor given by 

F,,(x-y) = -r2x,,(x-y) . (2.10) 

This result has two important features that motivate 
our formal considerations below. First, we anticipate 
that the pairing potential responsible for superconduc- 
tivity will be proportional to the magnetic susceptibility. 

Second, the interaction will be attractive if the spins form 
a relative spin- triplet, and repulsive if they form a spin- 
singlet.^^ 

2. Order parameter field theory 

The above considerations motivate transforming the 
fermionic field theory to an effective one in terms of 
an order parameter field for spin-triplet superconductiv- 
ity, and a magnetization field. Since we are considering 
an electronic mechanism for superconductivity, a (rela- 
tively) high Tc can be achieved only if electron-electron 
interaction effects are in some sense strong. This implies 
that for a consistent theory we need to also worry about 
strong-coupling effects as described by normal-state self 
energies, i.e., an Eliashberg-type theory is needed. Wc 
will accomplish this by introducing an additional field 
related to normal-state effects. 

Because we will eventually use a saddle-point, or free 
energy variational, technique we will have to choose a 
superconducting order parameter as an input. In the 
ferromagnetic phase it seems likely that the triplet su- 
perconducting phase with the lowest energy will be the 
one with the spins aligned in the direction of the magne- 
tization. In what follows we assume that this is the case. 
We plan to also investigate other superconducting order 
parameters in the future, see Section V below. 

In Appendix A we show how field-theoretic techniques 
can be used to derive standard strong coupling theory for 
conventional superconductors. Here we essentially repeat 
that procedure for the order parameter field mentioned 
above, with one modification. Because the Cooper in- 
teraction in the spin-triplet channel must be generated 
by magnetic fluctuations, we can no longer use a saddle- 
point approximation for all the degrees of freedom. Fur- 
ther, if the magnetic transition is of first order, this also 
manifests itself as a fliictuation effect in this description. 
Below we will see that it is sufficient to keep magnetic 
fluctuations to one-loop order in order to derive a theory 
for magnetic-fluctuation induced spin-triplet supercon- 
ductivity. Since the magnetic transition is already un- 
derstood, we will model its first order nature in a simple 
way. Generally, our guiding philosophy will be to keep 
fiuctuations to leading order if they lead to a qualita- 
tively new effect that is relevant for the superconducting 
transition temperature, but to neglect them otherwise. 

The first formal step is to introduce the magnetization 
field, M, by performing a Hubbard- Stratonovich trans- 
formation on the triplet interaction term, Eq. (2.2e). The 



other order parameter fields arc introduced via a delta- 
function constraint as in Appendix A, and we integrate 
out the fermions.^^' This allows us to write the partition 
function, in analogy to Eq. (A4), as 

J D[M, g, A] e-^[^'^'^l , (2.11) 

with the action A given by 

A[M, G,A] = ^ Jdx (tr [^oGix, x)]f + Tr (A^) 

+^Tv ln(Go 1 + ^/2^a • M - AT) 

-J dx M(a;) •M(a;) . (2.12a) 

Here we have constrained the field $, Eq. (2.5c), to a 
bosonic field Q by means of the Lagrange multiplier A, 
and the superscript T denotes a transposition operation. 
Tr is an operator that traces over all degrees of freedom, 
including the continuous space-time index, while tr de- 
notes a trace over discrete indices only. (70,7) is the 
matrix-valued four-vector defined in Eq. (2.6), and we 
have defined the inverse free electron Green operator. 

Go ' = -dr + 70 (VV2me + m) • (2.12b) 

The Eqs. (2.11) - (2.12) are formally exact. To proceed, 
we make a number of simplifications and approximations. 
First, we will assume a magnetic state, on the ferromag- 
netic (FM) side of the magnetic phase boundary, with 
an average magnetization m in the z-direction. We thus 
write 

M(a;) = m v'rt/2 z + 6M{x) . (2.13a) 

Motivated by the considerations at the beginning of the 
current subsection, we assume a non-unitary triplet su- 
perconducting order parameter {"ipiix) 'tpi{y))?'^ In terms 
of the matrix elements of Q, which are isomorphic to 
those of this order parameter and its hermitian ad- 
joint are 

2{gis{x,y)) = F+{x-y)^0 , 
2{gsi{x,y)) =F{x-y)^0 . (2.13b) 

We also define anomalous self energies, or gap functions, 

{Ai3{x,y)) = A+{x-y)j^0 , 
(A3i(x,y)) = A(x-y)^0 . (2.13c) 

The normal state order parameters and self energies are 
defined analogously to those in Appendix A, except that 
isotropy in spin space is lost in the ferromagnetic state, 

2{gn{x,y)) = -2{gss{y,x)) = G^ix - y) , 
2{g22{x,y)) = -2{g44y,x))=Gi{x-y) , (2.13d) 



4 



{An{x,y)) = -{A33{y,x)) = ^^{x-y) 
{A22ix,y)) = -{Kii{y,x)) = ^^(x- y) 



(2.13e) 



Note that we are assuming that the equihbrium state is 
spatially homogeneous. For our present discussion, which 
will focus on the superconducting phase transition, this 
assumption is legitimate. However, within the supercon- 
ducting phase it cannot be correct in general, as we dis- 
cuss in Section V. 

As already stressed, the crucial interactions leading 
to the superconducting state and the most important 
normal-state self-energy effects arc d\ie to magnetic fluc- 
tuations. To take these fluctuations into account, we in- 
sert Eq. (2.13a) into Eq. (2.12a), expand to order {8Mf, 
and then integrate out the magnetic fluctuations. The 
partition function and action can then be written 



,^[e,A] 



(2.14) 



with 

A[g,h\ 



Y j dx(tr [-iQg{x,x)\f 



f Tr (A^) 
^TrlnxX' 



(2.15a) 

Here we have defined a 4 x 4 matrix Green operator, 



Go 1 + rtm73 - A^ 



(2.15b) 



iXA^j'yX-y) =SijS{x-y) 

GA{y,x)'yiGA{x,y)-fj 



, (i,i = 1,2,3) , 
(2.15c) 



Notice that G and x both depend on A, as we explicitly 
indicate with our notation. 

The action as written in Eqs. (2.15) does not yet allow 
for a saddle-point solution that describes a superconduct- 
ing state, i.e. one with a broken U{1) gauge symmetry 
that corresponds to particle number conservation. This 
can be seen by taking the; variation with respect to, say, 
5i3, which yields A31 = 0, which characterizes a non- 
superconducting state. Indeed, such a solution is more 
than just a saddle-point solution, as can be seen by for- 
mally integrating out ^13. However, while always a so- 
lution of the field theory, such a state is not necessar- 
ily stable. To allow for superconducting states with a 
spontaneously broken gauge symmetry, we write A as its 
average value plus a fluctuation. 



A={A) + SA . 



(2.16) 



We then formally integrate over the 6A fluctuations to 
obtain 



/z>[51 



(2.17a) 



with 



A[Q, (A)] = m dx (tr Mx, x)]f + Tr ((A) g) 



1 



--Tr lnG/^\ - - Tr InX(A) - 



'(A) 2 



2T 

G7.\ 



+ ln [ D[SA] exp | Tr (gSA) -h ^ Tr In 



( A)+^A 
-1 
(A) 



-iTVln^^i^ 

2 - -1 



X 



(2.17b) 



(A) 



Notice that (A13) = A+, which couples to ^31 in the 
Tr ((A)(7) term, acts like an external fleld that breaks 
gauge invariance. The electron system with such an ex- 
ternal field defines an ensemble that allows for a nonzero 
superconducting order parameter, with (A) (and (g)) to 
be determined self-consistcntly.^^ 

The Green function G(a> {x — y) = G{x — y) will play 
a central role in our description of the superconduct- 
ing phase transition, so we explicitly give it here. It 
is convenient to express the matrix elements of G in 
terms of a normal-state Green function g. In frequency- 
wavenumber space, with a four- vector k = (a;„,k), we 
define 



and a normalized inverse magnetic susceptibility tensor, with 



2mc 



jjL — ad , 



(2.18a) 



(2.18b) 



where a = (t, i) = (+, — ) is the spin projection index, 
and 6 = F^m refiects the Stoner splitting of the Fermi 
surface in the ferromagnetic phase. We also define a func- 
tion 

D{k)=g:f\k)g:^\-k)+A{k)A+{k) . (2.19) 
Finally we note the symmetry relations. 



A(fc) = -A(-fc) 
A+(fc) = -A(-fc)' 



(2.20a) 
(2.20b) 



with the asterisk denoting a complex conjugate. The 
matrix inversion problem for G decouples into one 2x2 
problem and two 1x1 problems, and we obtain 



Gn{k)^g^\-k)/Dik) , 
Gi3{k) = A{k)/D{k) , 
G22{k)=giik) , 

Symmetry gives 



(2.21a) 
(2.21b) 
(2.21c) 
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G'33(fc) = -Gn(-fc) 

G44(fc) = -G22(-fc) 

G3i(fc) = -Gi3(A:)* 



(2.22a) 
(2.22b) 
(2.22c) 



All other Green's functions are identically equal to zero. 

Finally, the inverse magnetic susceptibility Xi^l) = 
is needed. Physically, x represents the normalized mag- 
netic susceptibility of a 'reference ensemble' of electrons 
that are subject to an external magnetic field equal to 5, 
and whose self energy operator is given by So-. General 
symmetry arguments show that the form 



X-'W = [5i^xl\k) + (1 - 5,3)x^;+(fc) 
+(5i2-^2i)x^'_(fc) , 



(2.23a) 



with the subscripts L and T denoting longitudinal and 
transverse components, respectively, and the transverse 
inverse susceptibility has both diagonal and off-diagonal 
pieces. Notice that the latter are nonzero only in a phase 
with a nonvanishing magnetization (or in the presence of 
an external magnetic field). For later reference we give 
the explicit expressions that result from our Gaussian 
approximation for the magnetic fluctuations. With Eqs. 
(2.15), (2.13), and (2.12b) we obtain 



Xl\k) = l + T, / Gn((7)Gn((7-fc) 

Jq 

+G22 (?)G22 (9 - fc) - G31 (g)Gi3 {q - fc)] , 

(2.23b) 

XT\{k) = 1 + T, [ Gii(g) [622(9 + k)+ 622(9 - k)] , 

Jq 

(2.23c) 

XtI_ {k) = iTt f Gil (q) [G22 (q + k) - G22{q - k) . 

Jq 

(2.23d) 

where / = T Y,n ^ '^^/i'^'^Y ■ More explicit expressions 
for these quantities will be discussed in Sec. Ill D below. 



III. SADDLE-POINT THEORY FOR 
SUPERCONDUCTIVITY DUE TO MAGNETIC 
FLUCTUATIONS 

A. Saddle-point solution 

The next step is to solve the field theory defined by the 
Eqs. (2.17) in a saddle-point approximation. Notice that 
there is only one fluctuating field left, namely Q. For fixed 
(A) we therefore have only one saddle-point condition. 



5A/5g\,^ = 



(3.1) 



and we will be looking for a translationally invariant 
saddle-point solution, Qsp{x,y) = Qspix — y) ~ {Q){x — y) 



(cf., however, the remarks in Sec. V below). Equation 
(3.1) yields an expression for the normal self energies, 
Eq. (2.13e), 

T,^ix-y) = -{dAn)A{x-y) 

-T,6{x - y)Y,GA^ - y) . (3-2a) 
^l{x-y) = -{SA22)A{x-y) 

-T,6{x-y)Y,GAx-y) , (3-2b) 

with Gff the normal Green functions, Eq. (2.13d). For 
the anomalous self energies, Eq. (2.13c), on obtains 



A+{x-y) 
A{x - y) 



-(5Ai3)A(a;- t/) 
-(5A3i)A(a;-?/) 



(3.2c) 
(3.2d) 



Here we have defined a 'A-ensemble' that is characterized 
by a statistical weight given by the integrand in the last 
term in Eq. (2.17b), with G replaced by w {Q), 



(...)a = j Dm ... 

where 

B[5K] = Tr (^spM) 
1 



j Dm 



,S[5A] 



(3.3a) 



■ Tr In 



^-1 



Tr In 



X(A)-|-5a/X(A) 



(3.3b) 



In order to make the Eqs. (3.2) more explicit, we now 
calculate {5h)\ in Gaussian approximation. Expanding 
the 'action' i?[(5A] to second order in (5A, and performing 
the Gaussian integral, we obtain 



2G- 



^GlG-1 



(3.4) 



Here we have neglected the contribution of the last term 

in B[5A], since it represents corrections to mean-field su- 
perconductivity due to magnetic fluctuations. In agree- 
ment with our general philosophy of keeping only lead- 
ing fluctuation effects, we drop these terms. Inserted in 
Eqs. (3.2), Eq. (3.4) provides a relation between the self 
energies and the saddle-point Green function. The free 
energy in saddle-point approximation, Fsp = — TlnZgp 
now reads 



sp 



2TT,Jdx(^G„{x,x)j 



-TTr ((A) ^sp) - ? Tr InG"! + | Tr Inx"^ 



-Tin j D[SA] 



spy 2 

B[5A] 



(3.5) 



and it depends on only one independent parameter, e.g., 
(A). 
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B. Minimization of the free energy, and 
strong-coupling equations 

As the last step in our formal development we now re- 
quire that the physical value of (A) minimizes the saddle- 
point free energy. The condition 



6F,p/5{A) = 



leads to 

{5A)a 





s 






S{A) 





(3.6) 



(3.7a) 



Here we have used Eq. (3.4), and a symbolic scalar prod- 
uct notation 





s 






S{A) 





^^^"^^^'^^^ ~ 



S{A) 

xxj^{y-x) . (3.7b) 



Notice that {SA)a, G~^, and (A) are still 4x4 matrices 
in spinor space (as well as functions of space-time), while 
X is a scalar in spinor space. Again in agreement with 
our general philosophy, we have neglected the contribu- 
tion of the last term in Eq. (3.5), which would lead to 
superconducting order-parameter fluctuation corrections 
to mean-field superconductivity. Equations (3.2) and Eq. 
(3.7a) now form a closed set of equations for the various 
self energies. Due to their matrix nature, they are com- 
plicated, and for our purposes we restrict ourselves to 
the linearized gap equations. That is, we evaluate Eqs. 
(3.7a) and (3.2) only to linear order in the gap parameter 
A. This procedure is sufficient for deriving the supercon- 
ducting phase boundary in both the paramagnetic and 
the ferromagnetic phases. We find 



A(A;) = 2ri 



J a 



q)A{q)\g^{q)\' 



(3.8a) 



S|(fc) = -r,n + 2rt [ [XL(fc-9)5T(9) 

Jq 

+2xT,+(fc-g)<7i(g)] , (3.8b) 
I]4(/c) = -r,n + 2rt / [xL{k-q)gi{q) 

Jq 

+2xT,+(fc-9)ffT(3)] , (3.8c) 

with 91 and gi given by Eqs. (2.18). A comparison with 
Eqs. (AlO) shows that Eqs. (3.8) are an obvious gener- 
alization of the conventional Eliashberg equations to the 
case of superconductivity mediated by spin fluctuations, 
with the magnetic susceptibility essentially replacing the 
phonon propagator, and the inequivalency of the two spin 
projections taken into account. 

The susceptibilities that appear in the Eqs. (3.8) are 
given explicitly by Eqs. (2.23). It is obvious from our 



derivation, however, that these explicit expressions are 
just a particular approximation for the exact suscepti- 
bilities, and that the latter will appear in the theory if 
one takes loop corrections into account. We will there- 
fore consider the x as physical susceptibilities that we 
model in the next subsections, using general physical ar- 
guments. We also note that in the paramagnetic phase, 
XL = XT. However, in the ferromagnetic phase the two 
susceptibilities are fundamentally different, as we will dis- 
cuss below. 

Finally, we have yet to discuss the magnetization as a 

function of the control parameters that drive the mag- 
netic transition. This is obtained by minimizing the free 
energy with respect to m, in order to find the physical 
values of the magnetization, 



6Fsp/5m = 



(3.9) 



A detailed theory of the itinerant ferromagnetic phase 
transition at low temperatures has been given elsewhere.^ 
In this paper, where we are interested in the supercon- 
ducting transition, we restrict ourselves to treating the 
magnetization in mean-field theory. That is, we neglect 
the fluctuation effects that are represented by the last 
two terms on the right-hand side of Eq. (3.5). We then 
obtain the magnetic mean-field equation of state in the 
form 



(3.10) 



We will discuss this expression explicitly in Sec. HID 2 
below. 



C. Gap equation for p-wave superconductivity 

So far the symmetry of the superconducting order pa- 
rameter has not been specified. To do so, we expand 
A(A;) in Legendre polynomials. 



A(fc) = A(iw„,|k|,cos0) =^A«(ia;„,|k|)Pi(cos^) , 

(3.11) 



where 9 is the azimuthal angle on the spherical Fermi 
surface. Since A(fc) = — A(— A;), the A^'^ are odd (even) 
functions of the frequency for even (odd) values of /. If 
we insert this expansion in Eq. (3.8a), we see that in the 
static limit only odd values of I contribute. Of the con- 
tributions with odd angular momenta, I = 1 will be the 
strongest for phase space reasons. We thus specialize 
to the p-wave case, and use the notation A^^^ = A. We 
further split off the frequency dependence of the normal 
self energy and write, as in Eliashberg theory,^^ 



Gn^^{(l)=i^nZ'({q)-Uf 



(3.12a) 



with 
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Z|(9) = l-[S|(ia;„,q)-I]t(0,q)]Aw„ , (3.12b) 



1. Paramagnetic phase 



and 



lq,T =<eq,T+5^T(0,q) 



(3.12c) 



The normal state Green functions in Eq. (3.8a) then pin 
the momentum to the up-spin Fermi surface, and we ob- 
tain for A^^\iuin,kp) = A{icOn) the equation 



(3.13a) 



where iuj„i = iujmZ-[{iu)rn)- In the same way we obtain 
an equation for Zf, 

w„ [1 - Z^{iuJn)\ = -TrT^sgn{wm) [dl{iuin - iiOm) 

+2d^{iu}n — iuJm)] ■ (3.13b) 
The kernels in these integral equations are given by 

dl{inn) = TtNl dxx{l-xy2)xj.{xkl,inn) , 
Jo 

(3.14a) 

dl{tnn)^TtNl dxxxi.{xkl,i^ln) , (3.14b) 
Jo 

d^{in„) = TtNl dxxxT{xkl,in„) . (3.14c) 

J X- 

Here iVp and are the density of states per spin at the 
Fermi level in the paramagnetic phase, and the density 
of states at the Fermi level for up-spin electrons in the 
ferromagnetic phase, respectively, and x± = 1 ± kp/kp. 



D. Explicit model for the longitudinal spin 
susceptibility 

For an explicit evaluation of our strong-coupling theory 
we need the longitudinal and transverse magnetic suscep- 
tibilities as input. For the determination of the super- 
conducting Tc, or the phase diagram, it suffices to know 
these susceptibilities in the normal-metal state. This is 
equivalent to the need of the phonon spectrum as input 
in ordinary Eliashberg theory. As in the latter case, one 
would ideally want to use experimental data as input, 
but this information is not available for the materials 
of interest. We therefore resort to theoretical modeling, 
with the aim of developing theoretical expressions that 
correctly reflect the relevant qualitative physical effects. 
Since the latter are different in the paramagnetic and fer- 
romagnetic phases, respectively, we will treat these two 
cases separately. 



In the paramagnetic phase, spin rotational invariance 
implies xl = XT = X- Equation (2.23b) gives an explicit 
expression which, after substituting the Green functions 
G and using Eqs. (3.8), results in an integral equation 
for xl- This is the result of a particular approximation 
we made in Sec. II, that integrated out the magnetic 
fluctuations in a Gaussian approximation. This partic- 
ular procedure results in a fairly elaborate treatment of 
the denominator of %, which takes into account the self 
energies of the Green functions, while the numerator is 
approximated by a constant, neglecting its wavenumber 
and frequency dependence. A slightly different Gaussian 
approximation was derived in Rcf. 30 by expanding about 
a Stoner saddle point within a theory focussing purely on 
magnetism. This approximation neglects the self energy 
S(j in the Green function g^, Eq. (2.18a), but keeps the 
wavenumber and frequency dependence of the numerator 
of X at the same level as that of the denominator. In the 
static, long-wavelength limit, and not too far from the 
magnetic transition, both of these zero-loop approxima- 
tions reduce to the Landau expression^^ 



X^?^(k,^0)=x^^^(k,^0) = 



(0)/ 



1-t 



t + 6(k/2fcF)2 



(3.15) 



Here t = 1 — 2A^Frt is the dimcnsionlcss distance from the 
magnetic critical point in mean-field approximation. The 
value of the coefficient b of course depends on the model. 
If one neglects the self energy in the Green function in 
Eq. (2.23b), or in the approximation of Ref. 30, b = 1/3, 
but in general there is no reason to prefer this value of b 
over any other value of order unity. 

We have used both the model from Ref. 30, and the 
Landau model, Eq. (3.15), to calculate the coupling con- 
stants dj^'rj,, Eqs. (3.14), and the resulting superconduct- 
ing Tc in a McMillan approximation, see Sec. IV below. 
We have fomid that for reasonable values of b the results 
are qualitatively the same. In what follows, we will use 
the Landau approximation, Eq. (3.15), with 6 as a free 
parameter. We emphasize that this approximation cor- 
rectly reflects the long-wavelength behavior of the mag- 
netic susceptibility in the paramagnetic phase. The free 
parameter b just reflects deviations from the free electron 
model that one would expect in any real material. 



2. Ferromagnetic phase, zero-loop order 

In the ferromagnetic phase, any reasonable zero-loop 
approximation for the susceptibilities yields again a Lan- 
dau expression for the longitudinal static susceptibility 
in the long-wavelength limit. We evaluate Eq. (2.23b), 
again neglecting the self energy, and find^^ 



x[?Hk,z0) = 



1-t 



aL|i|+&L(k/2/cF)2 



(3.16) 
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where we have introduced two free parameters, ol and 
6l- In ordinary Landau theory, ol = 2. However, if one 
keeps the electron number density fixed, as is the case ex- 
perimentally, rather than the chemical potential, then in 
general Ol ^ 2. From the point of view of Landau theory, 
this can be miderstood as follows. Since t = l — 2N(fj,)Tt, 
and the chemical potential fx depends on the magnetiza- 
tion, the coefficient of the magnetization squared in the 
Landau free energy depends itself on the magnetization, 
which causes the deviation of the pref actor of\t\ from the 
usual value of 2. 

The actual value of ol is model dependent, and de- 
pends in particular on the equation of state, since the 
latter determines the relation between |t| and the mag- 
netization. If we again neglect the self energy in the 
Green function on the right-hand side of Eq. (3.10), the 
latter yields 



3/2 



N3/2 



(3.17a) 



This needs to be supplemented by the condition that the 
electron number density rie is fixed, which is 



67rV = (2me)3/2 {fx + 6f^'' + {f, - 6) 



3/2 



(3.17b) 



This determines the chemical potential /i as a function of 
the magnetization. Equations (3.17) together determine 
the relation between S and \t\. 

From Equations (3.17) and Eq. (2.23b) one obtains 
ah = 1/2. Different Gaussian approximations for x yield 
different values. For instance, Ref. 30 in conjuction with 
Eqs. (3.17), and putting Fg = Ft, yields = 5/4, and 
so does the RPA expression used by Fay and Appel.^ In 
order to be able to compare with the latter's calculation 
of Tc, we adopt al = 5/4 in what follows. Other values 
of 0(1) do not lead to qualitatively different results. 

In the transverse channel, a Ward identity guarantees 
that X diverges in the limit of small wavenumbers and fre- 
quencies. This reflects the existence of magnetic Gold- 
stone modes or magnons. Our mean-field equation of 
state, Eq. (3.10), or Eq. (3.17a), is consistent with our 
Gaussian treatment of the magnetic fluctuations, Eqs. 
(2.23), in the sense that they form a conserving approx- 
imation that correctly reflects the Goldstone modes. In- 
verting the susceptibility tensor given by Eqs. (2.23), 
dropping the self energy in the Green functions, and us- 
ing Eq. (3.17a), we find for Xii(^) = X'T,+{k) in Gaussian 
approximation, 3 ^ 



)- ( 



/4eF + (5/2eF)6T(k/2fcF)2 
1 



inn/4eF - ((5/2eF)&T(k/2A:F)' 



(3.18a) 



and for Xi2{k) = xt,-(A:), 



XT^_(k,in„): 



i(5/4eF 



1 



-I- 



V-t V«^^n/4eF + (V2eF)&T(k/2fcF)2 
1 

iQ„/4eF - ((5/2eF)6T(k/2fcF)2 



(3.18b) 

Here we give the dynamical siisceptibility, since we will 
need it later. In our approximation, 6t = 1/3, but in 
general it is some parameter of order unity, and in general 
it will be different from 6l in the ferromagnetic phase. 



3. Ferromagnetic phase, one-loop order 

In contrast to the situation in the paramagnetic phase, 
in the ferromagnetic phase the zero-loop approximation 
for the longitudinal susceptibility, Eq. (3.16) is not qual- 
itatively correct. The reason for this is as follows. In 
a Heisenberg ferromagnet, the transverse spin waves or 
magnons couple to the longitudinal spin fluctuations, 
and hence contribute to the longitudinal susceptibility 
Xl-'^^ At nonzero temperature, and for dimensions d < 4, 
this coupling actually leads to a magnetic susceptibil- 
ity that diverges everywhere on the coexistence curve. 
More generally, and more importantly for our purposes, it 
causes the zero-loop or Landau expression for the longitu- 
dinal susceptibility, Eq. (3.16), to be qualitatively incor- 
rect in the ferromagnetic phase. Going to one-loop order, 
on the other hand, produces a functional form of xl that 
is exact at small wavenumbers. This one-loop contribu- 
tion can be calculated in a variety of ways. Suppose we 
expand the action, Eq. (2.12a), in powers of M, to quar- 
tic order. Then each of the four M flelds in the quartic 
vertex can represent either an average magnetization, or 
a magnetization fluctuation, see Eq. (2.13a). This leads 
to a one-loop contribution to the two-point M-vertex, or 
the inverse magnetic susceptibility, that has the structure 
shown in Fig. 2. Equivalently, an expansion in powers of 





FIG. 2. (a) A contribution to the quartic magnetization 
vertex, with the dashed line representing an average magne- 
tization, and the solid line a magnetization fluctuation, (b) 
A resulting contribution to the magnetic susceptibility. 

(5M leads to a cubic vertex, and to a diagram with the 
same structure as in Fig. 2(b). Alternatively, one can use 
a nonlinear sigma- model to derive the same contribution, 

as was done in Ref. 21. Either way one obtains for the 
one-loop contribution to the longitudinal susceptibility'^"'^ 

2Ft 



^2 



(0) 

Xl 



Jq 



X?V(^-9)xV7+(9) 



(0) 
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+ X?l(fc-«)X^°1(9)' 



x'^Hk) . (3.19) 



Here X^'± ^'"'^ zero-loop transverse susceptibility 

tensor elements, Eqs. (3.18), and Xi!"^ is the zero-loop 
longitudinal susceptibility, Eq. (3.16). The integral in 
Eq. (3.19) leads to XL(k,iO) oc \k\'^-^ for the leading 
wavenumber dependence of Xl in < 4, which is the ex- 
act leading behavior for small wave numbers."''^ We note 
that there is no analog of the diagram shown in Fig. 2(b) 
in the paramagnetic phase, as there is no cubic magneti- 
zation vertex. It is therefore reasonable to use the zero- 
loop expression for the longitudinal susceptibility in the 
paramagnetic phase, and the one-loop expression given 
by Eq. (3.19) in the ferromagnetic one. 

We now have explicit model expressions, in terms of 
integrals, for the magnetic siisc;eptibilitieH that we need 
as input for our theory of superconductivity. In the para- 
magnetic phase, the simple expression given in Eq. (3.15) 
suffices. In the ferromagnetic phase, we will use Eq. 
(3.18a) for the transverse susceptibility, and XL'''(k, iO) + 

x[)^(k,iO), Eqs. (3.16) and (3.19), with ol = 5/4, for the 
longitudinal one. This model contains two parameters of 
0(1), viz. 6l and 6t, and it reflects the exact functional 
form of the susceptibilities in the long- wavelength limit. 



IV. MCMILLAN-TYPE SOLUTION OF THE GAP 
EQUATION 

The strong coupling equations, Eqs. (3.13), have the 

same structure as Eliashberg theory for phonon-induced 
superconductivity. There are well-developed algorithms 
for solving these equations. It is well known that any 
attempt to quantitatively calculate the critical tempera- 
ture requires a numerical solution of the equations, using 
detailed information about the kernels as input. In the 
case of superconductivity induced by magnetic fluctua- 
tions, the latter information is not available, and we will 
have to use the schematic model for the susceptibility 
that was developed in the last section. We will therefore 
calculate the superconducting Tc by means of a simple 
expression of McMillan or Allen-Dynes type.'^^ Since we 
are interested in the structure of the phase diagram, and 
in particular in the relative values of the superconduct- 
ing Tc in the paramagnetic and ferromagnetic phases, 
respectively, these approximations will give the desired 
information provided we treat the two phases on equal 
footing. As we have explained in the last section, this is 
the case if we use the susceptibility model represented by 
Eqs. (3.15), (3.18a), and (3.19), respectively. For a rough 
order-of-magnitude estimate of the absolute value of Tc, 
see Sec. VA below. 



A. Tc formula 

McMillan's approximate solution of the linearized 
strong-coupling equations results in a critical tempera- 
ture 



Tc = To(i)e 



-{l+dl+2d°,)/dl 



(4.1) 



Here dj^ t ^^d are the zero-frequency values of the 
functions defined in Eqs. (3.14), and To(t) is a temper- 
ature scale that we will define and discuss below, d^, 
which involves an integral over a zero-loop susceptibility, 
is easily calculated analytically. For d!^^ we proceed as 
follows. The momentum integral of the convolution in 
Eq. (3.19) can be done analytically. For the external legs 
in that expression, we replace Eq. (3.16) by a step func- 
tion 0(|k| — 2k-pXc), with = \/5\t\/4h[. The momen- 
tum integration in Eq. (3.14b) can then also be performed 
analytically. This leaves us with the frequency summa- 
tion in Eq. (3.19), which needs to be done numerically. 
We obtain 



4 



Stt 



1 



52, ((5/eF)2 S+ S 



N 



7^2/0 + 2^/(0„,yc) 



71=1 



(4.2a) 



Here 6+ = y/l + S/ep, TJc = Min(l, Xc), and 



a-m + y 



+ ln2 



flm -y\ , / a„ 
arctan 



1 



arctan 



In 1 



am + y\ 
(am + yf 



In 1 



\am - y 

i2 



yf 



6(y - am)i^ 



y-a„ 



(4.2b) 



with Qm = \/^m/bTS. N in Eq. (4.2a) is given by TV = 
r2o/27rTc, with f2o a frequency cutoff that is on the order 
of 2&T^- Equation (3.14a) has an additional factor 1 — 
in the integrand, which we approximately take into 
account, in the same spirit as our approximation for the 
external legs, by writing 



1 



1 Tc 



6| (5/eF)2 S+ 6 



N 

2j2fi^n,yl) 



n=l 



(4.2c) 



with yl = Min(a;c, I/2V2). 

Finally, we need to specify the temperature scale To{t). 
Following Refs. 6, 35, we use the prefactor of |t| in Eqs. 
(3.15) and (3.16) as a rough measure of the magnetic 
excitation energy. 



To{t)=To [Q{t)t + e{-t)5\t\/4], 



(4.3) 



with To a microscopic temperature scale that is related 
to the Fermi temperature, for a free electron model, or 
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a band width, for a band electron model. This choice of 
To{t) qualitatively reflects the suppression of the super- 
conducting Tc near the ferromagnetic transition due to 
effective mass effects. 

Equations (4.1) - (4.3) now form a closed set of tran- 
scendental equations for the superconducting Tc, which 
we have solved by numerical iteration. We discuss the 
results in the next subsection. 



B. Results 

We have solved the Tc equations (4.1) - (4.3) for 
various values of the parameters, and with various 
implementations of our approximations. Specifically, 
we have also used RPA expressions for the magnetic 
susceptibilities^^'^° instead of the simple Landau expres- 
sions given above. While this makes the calculations 
substantially more complicated, the results were qualita- 
tively the same. We therefore only show results obtained 
by using the Landau model. We have also modified the 
upper frequency cutoff, and again have found only quan- 
titative effects. Since Tc depends exponentially on the 
coupling constants, the quantitative effect of changing 
the cutoff is substantial. However, the relative values 
of Tc in the paramagnetic and ferromagnetic phases, re- 
spectively, depend only very weakly on the cutoff. We 
therefore use flo = 26t^ throughout. Finally, in the fer- 
romagnetic phase it is illustrative to relate the supercon- 
ducting Tc to the magnetization in addition to t, so we 
need again the magnetic equation of state. Within our 
mean-field model for the magnetization, Eqs. (3.17) it is 
useful to introduce a parameter r] via^^ 



i = l-(l + 3r;2)i/3/(i + ^2/3) 



(4.4a) 



The magnetization, in units of fiB^ with the Bohr 
magneton and n the electron number density, is related 
to T] by^^ 



m/flBTle = 3??(1 + ??V3)/(1 + 377^ 



(4.4b) 



Finally, we take into account theoretical^ and 
experimental^ evidence that the quantum ferromagnetic 
transition in clean itinerant electron systems is of first 
order at least in some systems. This is because fluctua- 
tion effects first appear at one-loop order lead in general 
to a negative term in the Landau free energy, i.e., one 
has a fluctuation-induced first order transition.^'^ A very 
simple way to take this into account in the present con- 
text is to just impose a minimum value imin on \t\, which 
corresponds to a minimum value of the magnetization in 
the ferromagnetic phase. 

Figure 3 shows the superconducting Tc as a function 
of t with t,„in ~ 2.66 X 10^'*, which corresponds to a 
discontinuity in the magnetization equal to 6% of the 
saturation value, 6 = 6l = 0.23, and 6t = 0.25. Also 
shown is the magnetization m in the ferromagnetic phase 



in units of the saturation magnetization /iefio- Tc is mea- 
sured in units of the characteristic temperature Tq that 
is given by either the Fermi temperature or the band 
width, depending on the model considered. The solid 
curves show the result in the paramagnetic and ferro- 
magnetic phases, respectively, with the former scaled by 
a factor of 133.3 (right-hand scale) compared to the latter 
(left-hand scale). The dotted curve in the ferromagnetic 
phase (also scaled by a factor of 133.3, right-hand scale) 
represents the result that is obtained upon neglecting the 
mode- mode coupling effect. The zero-loop results in the 
paramagnetic and ferromagnetic phases are very close to 
one another, and also very close to the result of Fay and 
Appel.^ Upon inclusion of the mode-mode coupling ef- 
fect at one- loop order, however, the maximum value of 
Tc in the ferromagnetic phase is enhanced by a factor 
of more than 100. As we have explained above, there is 
no analogous contribution in the paramagnetic phase, so 
this comparison is physically sensible. We also reiterate 
that one should not take the absolute Tc values very seri- 
ously. However, the relative comparison we expect to be 
reliable. We thus find a pronounced asymmetry between 
the paramagnetic and ferromagnetic phases. In the case 
of UGe2, where the observed masimum Tc in the ferro- 
magnetic phase is about 500 mK, we thus predict that 
in the paramagnetic phase no superconductivity should 
be expected above temperatures of at most about 5 mK. 
This is in agreement with the experimentaUy observed 
absence of superconductivity in the paramagnetic phase. 
The second maximum in the Gaussian Tc in the FM phase 
we will discuss below. Fig. 4 shows the result with the 
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FIG. 3. Superconducting Tc as a function of the distance 
from the critical point t, and the magnetization m. The solide 
curve in the PM phase and the dotted curve in the FM phase 
(right scale) show the zero-loop Tc scaled by a factor of 133.3, 
while the solid curve in the FM phase (left scale) represents 
the one-loop result. See the text for parameter values and 
further explanation. 

same parameters, except that tmin = 6.5 x 10~^, corre- 
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spending to a magnetization discontinuity equal to 29% 
of the saturation value. Figures 3 and 4 correspond to 
systems with a magnetic phase transition that is weakly 
and strongly first order, respectively. 
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FIG. 5. Same as Fig. 3, but for different parameter values. 
See the text for further explanation. 



FIG. 4. Same as Fig. 3, but with a different value of tmin, 
see the text. 

We now discuss the effect of varying the parameters 6l 
and br- With 6 = 6l = &t = 0.5 and t^in = 1.7 x 10"^ 
we obtain the result shown in Fig. 5. 6 = 6l = &t = 1-0 
and tjnin = yields Fig. 6. The left and right-hand scales 
differ by a factor of 50 in both figures. While the (un- 
physical) zero-loop result in the ferromagnetic phase is 
very sensitive to the parameter values, we see that the 
enhancement of the (physical) one-loop result over the 
in the paramagnetic phase is rather robust. However, the 
position of the maximum of Tc changes compared to Fig. 
3; in Fig. 6 it occurs at the point where the magnetization 
reaches its saturation value, and in Fig. 5 there are two 
maxima of about equal hight. The reason is as follows. 
As one approaches the magnetization saturation point 
from low magnetization values, the transverse coupling 
constant c?^ vanishes, and remains zero in the saturated 
region. Effectively, the Heisenberg system turns into the 
Ising model discussed in Ref. 19. If the longitudinal cou- 
pling constant dj^ still has a substantial value at that 
point, then this leads to an increase in Tc. This is a very 
strong effect in the zero- loop contribution, see Figs. 5 and 
6, and the effect qualitatively survives in the one-loop re- 
sult. If, however, is already very small, then dip going 
to zero has only a very small effect on Tc, as is the case in 
Figs. 3 and 4, although the effect is still visible in the zero- 
loop result. Which of these two cases is realized depends 
on the parameter values. This observation provides an 
explanation for the difference in the observed phase dia- 
grams in UGe2 and ZrZn2, respectively. UGe2 reaches its 
saturation magnetization relatively close to the magnetic 
phase transition, at a pressure well above ambient pres- 
sure. ZrZn2, on the other hand, has not nearly reached 
its saturation magnetization at ambient pressure. The 
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FIG. 6. Same as Fig. 3, but for different parameter values. 
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above discussion suggests that, as a consequence, the 
superconducting Tc in ZrZn2 is likely to monotonically 
increase from its onset as one goes deeper into the fer- 
romagnetic phase by decreasing the pressure, while in 
UGe2 one expects the superconducting Tc to peak close 
to the magnetic phase boundary, leading effectively only 
to a pocket of superconductivity, as shown schematically 
in Fig. 1(a). This is indeed what is observed in the ex- 
periments, and the present model provides a qualitative 
explanation. Consistent with the experimental observa- 
tion that the magnetic transition in ZrZn2 is at most very 
weakly first order, we have used tmin = in Fig. 6. 
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V. DISCUSSION 

We conclude with a summary of our results, and then 
briefly discuss several open questions. 

A. Sumrticiry, and additional remarks 

We have made two contributions with this paper. The 

first one is general in nature, and consists in a field- 
theoretic framework for deriving strong-coupling equa- 
tions for superconductivity, regardless of the mechanism 
that causes superconductivity. This framework is a nat- 
ural extension and generalization of Landau's theory of 
phase transitions, and it has several advantages, both 
philosophical and practical ones, over the many-body di- 
agrammatic formalism that is usually used to describe 
superconductivity. One advantage is the fact that the 
field theory allows for introducing several order param- 
eters simultaneously in a very natural and transparent 
way. In the present case wc have used this to obtain a 
theory of the interplay between ferromagnetism and spin- 
triplet, p-wave superconductivity that deals with both or- 
dering phenomena on equal footing, and within one uni- 
fied framework. No input from other theories has been 
necessary. We have restricted our description of the su- 
perconductivity to the mean-field level, but it is obvious 
how to include fluctuations within this framework. 

In this context wc note that for the pairing mechanism 
we have assumed, or for any purely electronic mecha- 
nism, there is no separation of energy scales, and hence 
no analog of Migdal's theorem. As a result, Eliashberg 
theory cannot be expected to produce quantitative re- 
sults, in constrast to the case of phonon-induced super- 
conductivity. A rough order-of- magnitude estimate of 
may nevertheless be of interest. Fay and Appel^ have 
given detailed estimates of the zero- loop Tc for ZrZn2. 
These ranged from 100 /uK to 500 mK, depending on 
parameter values that are not well known. With a mid- 
range value of 5 mK, which is consistent with the ob- 
served absence of superconductivity in the paramagnetic 
phase, our enhancement mechanism yields a maximum 
Tc in the ferromagnetic phase of about 500 mK, which 
is the same order of magnitude as observed experimen- 
tally. Regardless of how seriously one is willing to take 
such estimates, it is important to note that any correc- 
tions due to the inapplicability of Migdal's theorem are 
expected to equally affect both the paramagnetic and the 
ferromagnetic phases. They therefore will not spoil the 
relative comparison of the two phases that has been our 
main objective. 

We have also restricted ourselves to discussing the lin- 
earized gap equation, but our formalism is more gen- 
eral, and includes the feedback of a nonvanishing super- 
conducting order parameter on the magnetic properties. 
Since the magnetic properties arc responsible for the su- 
perconductivity in the first place, this feedback will have 



to be included self-consistently in a theory of the super- 
conducting phase. These points will be pursued in future 
publications. We further note that we have assumed the 
superconducting transition to be of second order. First 
order scenarios are possible, see Ref. 38. 

Our second contribution is specific, and consists of 
a calculation of the phase diagram for p-wave, spin- 
triplet superconductivity that is driven by ferromagnetic 
spin fluctuations and coexists with ferromagnetism. The 
most striking qualitative feature of the result is the pro- 
nounced asymmetry of the superconducting critical tem- 
perature in the ferromagnetic and paramagnetic phases, 
respectively, which results from the contribution of the 
magnons in the ferromagnetic phase to the pairing mech- 
anism. The results are in qualitative agreement with the 
experimental observations in UGe2 and ZrZn2. 

B. Outlook 

The most intriguing open questions concern the na- 
ture of the magnetic-fluctuation induced superconduct- 
ing state, and an understanding of the phase diagram for 
all temperatures, magnetizations, and magnetic fields. 
One interesting question is whether there is more than 
one superconducting phase as a function of temperature, 
magnetization, and external magnetic fields. Since the 
pairing mechanism is expected to be electronic in ori- 
gin, and itself sensitive to superconductivity, it is easy 
to imagine additional superconducting states appearing 
inside the superconducting phase as the temperature is 
lowered. In particular, the Cooper pairs that form in 
the superconducting state produce an additional inter- 
nal magnetic field. If this field acts like an external field 
for the electrons, then the fermionic spin waves, i.e., the 
ferromagnetic Goldstone modes, will acquire a mass from 
this effect, leading to a smaller fluctuation contribution 
to the Tc mechanism. This suggest another phase in- 
side the superconducting one, where Cooper pairs anti- 
aligned to the magnetization might also exist. In general, 
we also note that the concepts of transverse and longitu- 
dinal critical magnetic fields needs to be worked out for 
these superconducting states. 

To understand the precise nature of the superconduct- 
ing state, another fundamental problem must be ad- 
dressed. In the theory presented in Section II, the only 
coupling of the internal magnetic field in the ferromag- 
netic phase to the superconducting order parameter is 
through the Zeeman-like term that leads to the Stoner 
splitting of the Fermi surface, and through the Goldstone 
modes (which in some sense are also due to the Zecman 
term) that exist in the ferromagnetic phase. An obvi- 
ous question is then whether there arc orbital effects in 
the Cooper channel due to the internal magnetic field. If 
there are such effects, they will lead to an inhomogeneous 
superconducting state, viz., to some type of Abrikosov 
flux-lattice state. This magnetic field effect is of relativis- 
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tic origin and has been neglected in the current theory. 
It has not yet been discussed theoretically for pairing 
mechanisms that are of electronic origin and sensitive to 
internal magnetic field effects. 

Finally, our theory of the superconducting phase tran- 
sition is a mean-field treatment. The validity of this type 
of approximation in estimating a superconducting criti- 
cal temperature is not obvious. We note, however, that 
one would expect any fluctuation-induced degradation of 
the superconducting to be larger in the paramagnetic 
phase than in the ferromagnetic one, since for triplet su- 
perconductivity, the direction of the internal magnetic 
field caused by the magnetization will lead to an addi- 
tional phase coherence of the triplet Cooper pairs in the 
magnetic phase. The relative importance of this effect 
on the phase diagram needs to be investigated. 
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APPENDIX A: CONVENTIONAL ELIASHBERG 
THEORY 

In this appendix we show how to derive conventional 
Eliashberg theory by means of field-theoretic methods. 
This is an extension to strong-coupling theory of previous 
field-theoretic derivations of BCS theory.^^ 

We start with the fermionic field theory defined in Sec. 
II A, but add harmonic phonons, 

Sph = j dxdy (p{x)D^^{x-y)f{y) , (Ala) 

and an electron-phonon interaction, 

Se-ph = g dx(p{x)^-4)aix)llJa{x) ■ (Alb) 
a- 

Here D^^ is the phonon vertex function whose inverse, 
D{x — y) = D{y — x), is the phonon propagator, and g 
is the electron-phonon coupling constant. For simplicity, 
we consider a scalar phonon field (fiix); the phonon po- 
larizations can be restored in an obvious way if desired. 

We now proceed to integrate out the phonons. This 
produces an effective electron-electron interaction, which 
we decompose into particle-hole and particle-particle 
channels, and spin-singlet and spin-triplet contributions 
in each channel, as usual. In the particle-particle chan- 
nel, only the spin-singlet is nonzero, and in the particle- 
hole channel we neglect the spin-triplet contribution 



Sf ^, Eq. (2.2e). The electronic part of the action then 

reads 



S = - jdxY^ i){x) dr i){x) + Sq + 5p-^ 



1 cP-h I CP-P 
"'"'^s.e-ph '^s.e-ph 



(A2a) 



with 

c'P-r 

s,e— ph 



5f „ P„v, = I dx dy V'i (x) ipi (y) D(x - y) ipi (y) ipi (x) 



cP-ii _ ~y 

s,e— ph ^ 



■ dxdyY^ 



i'aix) il)„{y)D{x - y) 



X'^a'{y)i^a'{x) 



(A2b) 



(A2c) 



and So and ^P"^ from Sec. II A. 

Now we introduce the Nambu spinors 



{ i'lix) 



and the corresponding composite variables 
<^ij{x,y)=^i{x)^j{y) . 



(A3a) 



(A3b) 



The field $ is bilinear in the fermion fields, and hence 
its components commute with each other as well as with 
all other objects. $ is therefore isomorphic to a classical 
field G, and we can transform to a description in terms 
of these variables by exactly rewriting the partition func- 
tion, Eq. (2.1), as 



= J D[ij,i^] e^I'^-'^l J D[g] 6[g-^] 

= y e^[*'*l J D[g,A] e^[A(e-*)l 

= J D[g,A] e-^[^'^] . (A4) 



Here A is an auxiliary field that serves as a Lagrange 
multiplier, and we have defined an effective action 

A[g, A] = ^Tr ln(Go 1 - A^) + Tr (A^) 

+g'^ j dx dy ^12 {x,y)D{x- y) ^21 {y, x) 

dxdy [gii{x,y) D{x - y)gii{y,x) 

+g22{x,y) D{x - y) g22{y,x) 
-2gii{x,y) D{x - y) ^^22(2/, a;)] 

\2 



+ 



~^ J dx (tr(T3 C/(a;,a;))^ . (A5a) 



Tr denotes a trace over both discrete and continuous 
labels, while tr traces over discrete labels only. Gq^ is 
an inverse free-particle Green operator, 
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Go' = -dr + 



(^3 



2mp 



(A5b) 



with fT3 the third Pauh matrix. The matrix elements of 
Go are 



Go{x ~y) = 



_ fGo{x-y) 









-Go{y-x) 



, (A5c) 



with 



Go (a; -y) = ipcr{y))so 

= {x\ -dr + VV2me + n\y) , (A5d) 

the usual free particle Grcoii function. 

The standard strong-coupling, or Eliashberg, theory 
of superconductivity is just the saddle-point solution of 
the field theory given by Eqs. (A5). The saddle-point 
condition is 



^ SA ^ 



(A6) 



We are looking for homogeneous saddle-point solutions 
Gsp{x,y) =Q{x-y) and Asp(x,?/) = A(a; - y), and it is 
easy to see that the diagonal components of both Q and 
A are related. To express this fact, and to conform with 
standard notation, we write 



r( _ \-( G{x-y) -F+{x-y) 
y)-\_F{x-y) -G{y-x) 



K{x - y) 



E{x-y) A+{x-y) 
A{x-y) -^{y-x) 



(A7a) 
(A7b) 



As we will see, G and F , are the usual normal and 
anomalous Green functions, and S and A, A"*" are the 
normal and anomalous self energies, respectively. 
The first equality in Eq. (A6) gives 

S(a; -y)= [g^ D{x - y) - 2V, 5{x - y)] G{x - y) , 

(A8a) 

A+{x-y)^g^D{x-y)F+{x-y) , (A8b) 
A(x - y) = D{x - y) F{x - y) . (A8c) 

The second equality yields 

G{x -y) = Go(x - y) 

+ j dx' dy' Go{x - x') i:{x' - y') G{y' - y) 

+ j dx' dy' Go{x - x') A{x' - y') F+{y' - y) , 

(A9a) 

F+{x -y)= j dx' dy' Go(x' - x) E(y' - x') F+{y' - y) 

+ J dx' dy' Go{x' - x) A+{x' - y') G{y' - y) , 

(A9b) 



F{x -y)= j dx' dy' Go{x - x') E{x' - y') F{y' - y) 

+ j dx' dy' Go{x - x') A{x' - y') G{y - y') . 

(A9c) 

The Eqs. (AS) and (A9) are the standard EHashberg 
equations for conventional, spin-singlet, phonon-induced 
superconductivity.^^ Vg plays the role of the Coulomb 
pseudopotential.''" 

The linearized gap equation, which determines the su- 
perconducting transition temperature, is obtained by ex- 
panding the Eliashberg equations to linear order in A. 
In Fourier space, with frequency-momentum four- vectors 
k = (ia;„,k), one finds for the linearized gap equation 

A{k) =g^ [ D{k - q) G{q) A(g) G{-q) , (AlOa) 

Jq 

with a normal Green function 

G(fc) = l/[Go^(fc)-S(fc)] , (AlOb) 
in terms of a normal self energy 

E(/e) =g^ f D{k - q) G{q) - V^n^ , (AlOc) 

Jq 

Here Jg = Tj2nI dci/{2T^ f, and = 2TrG is the elec- 
tron number density. 
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